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Abstract 

We examine the effects of relic dark matter annihilations on big bang nucleosynthesis (BBN) . 
The magnitude of these effects scale simply with the dark matter mass and annihilation cross- 
section, which we derive. Estimates based on these scaling behaviors indicate that BBN severely 
constrains hadronic and radiative dark matter annihilation channels in the previously unconsid- 
ered dark matter mass region MeV < m x < 10 GeV. Interestingly, we find that BBN constraints 
on hadronic annihilation channels are competitive with similar bounds derived from the cosmic 
microwave background. 



1 Introduction 



Dark matter populates our universe, representing more than 80% of the total matter content. 
While astrophysical and cosmological measurements have elucidated the amount of dark matter 
(DM) in our universe as well as its role in galactic formation and distributions, we know very little 
about the nature of DM. For example, one of the few properties of DM that we can say anything 
with certainty is that dark matter has mass. But even the allowed mass values can span over 80 
orders of magnitude, 10 -22 eV < m x < 10 59 eV, with the lower bound coming from requirement 
of localization on galactic scales [1] and the upper bound from micro-lensing searches for compact 
halo objects [2]. 

Despite our ignorance, we remain hopeful that theoretical models can provide a framework for 
experimentally probing dark matter directly and indirectly. For example, in this work we consider 
dark matter to be a weakly interacting massive particle that started its cosmic history in thermal 
equilibrium and froze-out as the universe cooled. The assumption of thermalization provides two 
powerful statements on DM: first, it restricts the allowed mass region to 11 orders of magnitude, 
keV < m x < 100 TeV, with the lower and upper bounds coming from the requirement that dark 
matter be cold [3] and its annihilation unitary [3], respectively. Second, if DM is a thermal relic, 
its annihilations must freeze-out to reach the observed current abundance, providing an estimate 
on the strength of DM interactions. In particular, assuming annihilation is s-wave dominated, in 
order to meet the observed abundance, freeze-out requires a weak scale cross-section that is nearly 
independent of mass, (crv}th — 3 x 10~ 26 cm 3 /s. 

After freeze-out, relic annihilations that occur may have observational consequences in two 
different scenarios. First, present day annihilations may be indirectly observed by searching the sky 
for their annihilation products (for a review, see [5]). In the second scenario, relic annihilations can 
inject hadronic and/or electromagnetic energy that may alter events in our cosmic history, namely 
big bang nucleosynthesis (BBN) and recombination (CMB). Because the physics of nucleosynthesis 
and recombination are well understood, they offer particularly clean environments in which to 
probe non-standard effects, such as DM annihilation. In addition, changes to BBN and CMB by 
non-standard processes depend only on the type and rate of energy injection into the bath, and 
therefore leave BBN and CMB essentially decoupled from the high-energy details allowing relatively 
model-independent statements to be made. 

In this work, we focus specifically on how relic annihilations from thermal dark matter affect 
nucleosynthesis. However, it is useful to describe in general how injected hadronic and electromag- 
netic energy alter nucleosynthesis and recombination. During nucleosynthesis, injection of hadronic 
and/or electromagnetic energy alters the abundances of nuclei via (i) energetic nucleons and pho- 
tons dissociating nuclei and (ii) pions inter-converting background nucleons, thus raising the n/p 
ratio and therefore the primordial 4 He mass fraction. As for recombination, injected electromag- 
netic energy ionizes hydrogen atoms and therefore broadens the surface of last scattering. This 
results in scale dependent changes to the temperature and polarization power spectra, especially 
in the low I modes. 

Now that we live in an era of precision CMB and BBN measurements, quantifying the effect of 
energy injection has led to powerful constraints on non-standard effects such as decay of long-lived 
particles P El E], El HUl HH [121 H31 E] and dark matter annihilation P [ISl [TO1 [IS1 [J71 [IHl [111 EDI EIJ 
[22 l 123 1 [2^ 1 125] . in the case of dark matter annihilation, the primary result obtained from analysis of 
the injected energy on nucleosynthesis and recombination is an upper bound on the DM annihilation 
cross-section, (crv), as a function of the DM mass. In particular, detailed analysis of the CMB has 
constrained s-wave annihilation to hadronic or electromagnetic channels to be less than the thermal 
cross-section required for freeze-out, (av) < (crv)th, for m x < 10 GeV [201 E21 [Ml [25] . Studies of 
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s-wave annihilation on BBN have been performed by Hisano et. al. for m x > 100 GeV |16t [T7] and 
separately by Jedamzik and Pospelov for m x > several GeV |18j . while lower mass regions have 
remained unconsidered. 

From previous works |16l \T7\ [T8] that consider the effect of DM annihilation on BBN, the 
constraint on the annihilation cross-section appears to have a simple power law dependence on the 
DM mass, (av) oc m x , where the power for hadronic (electromagnetic) energy injection is 5 ~ 3/2 
(= 1). In this work, we explain these scaling behaviors and the range of DM masses for which they 
hold. Then we extrapolate the bounds to previously unconsidered DM masses and find that BBN 
constrains (av) < (av) t h for few GeV < m x < 10 GeV (30 MeV < m x < 75 MeV) for the case of 
hadronic (electromagnetic) energy injection. 

The structure of this paper is as follows: in section II we discuss how injection of hadronic and 
electromagnetic energy alters BBN, and in particular, how these effects scale with the dark matter 
mass and annihilation cross-section. In section III, we use these scaling laws along with results 
from precise numerical treatment [16] to estimate constraints on DM annihilation in a previously 
unconsidered low mass region. Our estimates indicate that BBN places a constraint on hadronic 
annihilation channels that is competitive and independent of the CMB constraint. We also discuss 
how these bounds change for annihilation to other Standard Model particles. Finally, in section IV 
we conclude and discuss future prospects for this work. 



2 Energy injection into BBN 

Injected energy from dark matter annihilations occurring after T ~ 1 MeV can alter the primordial 
abundances of the light elements. Details of the evolution of injected energy and its effects on 
nucleosynthesis are thoroughly discussed in [3J. Quantifying the change on primordial abundances 
boils down to calculating the spectrum of produced hadron/nucleus Hi. This spectrum, fjg t = 
driHil 'dEjii-, is a complicated function that depends on the dynamics of the injected energy, the 
nuclear network leading to the production/destruction of Hi, and the rate of DM annihilation. 
In this work, we focus on the dependence of this spectrum on the DM's mass and annihilation 
cross-section. Since production/destruction of Hi begins with an annihilation event, dfjji/dt is 
proportional to the DM annihilation rate, 

dfHi dn x 2 

Here, (av) is the thermally- averaged annihilation cross-section which we take to be s-wave dom- 
inated and time independent and n x is the number density of DM. For the range of DM masses 
considered, BBN occurs after freeze-out and therefore the number density of DM is fixed by obser- 
vational abundance 

= i K J_. (2) 

m x a(t)° m x 

Hence, the DM injection rate scales as 



dn x (av) 

IT oc 2 

at m z 



x 



This has a simple physical interpretation: as dark matter mass decreases, in order to fit the observed 
abundance there needs to be more DM particles, which in turn leads to more annihilation events. 
In addition to the DM annihilation rate, the spectrum of produced/destroyed Hi will depend on 
the DM mass through the dynamics of the type of energy injected, which we now consider. 
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2.1 Injection of hadronic energy 

Hadrons (p, n, 7r ± , etc.) resultant from DM annihilation can primarily alter BBN through proton- 
neutron interconversion (ir ) and hadro-dissociation (p, n) IS]- Interconversion of p and n by 
energetic pions before the formation of D and 4 He (T ~ 100 keV) can increase the n/p ratio set by 
v decoupling at T ~ MeV. Since almost all neutrons end up in 4 He (see, e.g. [26J), the increase of 
n/p leads to a larger 4 He primordial mass fraction, Y p . 

Energetic nucleons produced by DM annihilation alter primordial abundances through collisions 
with nuclei. The dominant effect of energetic neutrons and protons is to hadro-disassociate 4 He 
which in turn leads to an increased production of D, 3 H, 3 He, and 6 Li. Electromagnetic interactions 
with background photons and electrons tend to thermalize the nucleons, so hadro-dissociation 
only becomes efficient when the universe dilutes enough so that the electromagnetic interactions 
lose their stopping power. This occurs below T ~ 100 keV for neutrons and T ~ 10 keV for 
protons [9j. Because the primordial abundance of 4 He is much larger than other nuclides (for 
example, no/^jje ~ O(10 -3 )), BBN is more sensitive to enhanced production of other nuclides 
than to processes which alter the primordial 4 He abundance. Therefore, hadro-dissociation of 4 He 
leading to increased production of other nuclides — namely D, 3 H, 3 He, and 6 Li — provides the most 
stringent constraints from BBN. 

Since we wish to understand the change in a given nuclide's primordial abundance as a function 
of {av) and m x , we consider how the amount of hadro-dissociation depends on these parameters. As 
in equation (pQ), the amount of hadro-dissociation is proportional to the dark matter annihilation 
rate, dn x /dt oc (avj/nit. Since hadro-dissociation proceeds through energetic nucleons, further 
dependence on m x essentially comes from the number of energetic nucleons produced in an anni- 
hilation event, i.e. the nucleon multiplicity [BJ[7]0 Assuming annihilation into quarks, XX ~* QQi 
the nucleon multiplicity can be approximated by the multiplicity in e + e~ collisions at y/s = 2m x . 
To leading order in QCD, the average charged particle multiplicity has the following dependence 
on V« [ZZIEEIEH]: 

(rich) = aexp ( W^-Wlog (^) ) +c (4) 



where /3q = 11 — 2nj/3, A is the renormalization scale, and a and c are constants. Apart from overall 
normalization, for \fs = 2m x above several GeV, the above is well approximated by (\/s)°' 5zize ~ 



m°- 5±e with e ~ 0.05. 

A 



The dependence of the amount of hadro-dissociation on the rate of DM annihilation and the 
nucleon multiplicity implies that the Boltzmann term for production / destruction of nuclei Hi ap- 
proximately scales with DM as: 



dn Hi 
dt 



K m X m 2 3/2 W 

hadronic " L X TTly 



This scaling is expected to hold down to m x ~ few GeV, below which hadro-dissociation quickly 
goes to zero since an annihilation event can no longer produce nucleons. This behavior implies 
that constraints on (av) from BBN approximately scale as m x , which is confirmed by precise 
numerical calculations for m x > 10 GeV [HI El El] an d shown in figure [TJ 



technically, it is the spectrum of nucleons per annihilation event, dN^ /dEx, that is of interest. The dependence 
of this spectrum on y^i = 2m x comes from the low x = En/^/s behavior of the parton fragmentation functions 
Di(x, y/s). The dependence of the multiplicity on y/s shown in the text is inherited from dNN /cLEat's dependence 
on y/s [23G8]. 
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2.2 Injection of electromagnetic energy 

Electromagnetic showering of energetic photons and leptons in the early universe produces many 
photons [30j which can photo-dissociate nuclei and alter primordial abundances [U [9] . 

Photo-dissociation is a significant effect only when the energetic photons do not thermalize. 
Precise calculations of the evolution of electromagnetic cascades in the early universe can be found 
in the literature |30j . but the dominant effects are easily understood. Because of the small baryon-to- 
photon ratio, 77 ~ 0(1CT 10 ), photons will pair produce off background photons (7 + 7bk g — > e + + e ) 
if their energy is above the threshold energy, E t ^ ~ m^,/22T. Therefore, photo-dissociation only 
becomes efficient when the binding energy drops below the threshold for pair production. 

Of the possible photo-dissociation channels, destruction of 4 He dominates since He is the 
most abundant of the light nuclei by a few orders of magnitude. Since E b 4jj e ~ 20 MeV, photo- 
dissociation of 4 He does not become relevant until T < 0.5 keV. The dominant photo-dissociation 
channels of 4 He are to 3 H or 3 He plus a nucleon [8j. Therefore, photo-dissociation of 4 He most 
noticeably leads to overproduction of 3 He ( 3 H decays to 3 He with a half-life of about 14 years). 
Overproduction of D and 6 Li may also result from 4 He photo-destruction. In the case of deuterium, 
it can be directly produced in photo-spallation of 4 He [8 J or in capture of the neutron from 7+abkg — > 
3 He+n on a background proton. Non-thermal production of 6 Li may occur if the spallation product 
3 H or 3 He is captured on a background 4 He, 3 H + abk g — > 6 Li + n or 3 He + abkg — > 6 Li + £>- Because 
the observed abundance of 6 Li is very small, non-thermal production of 6 Li via photo-dissociation 
of 4 He can have a noticeable impact [31|, [32] . 

Assuming annihilation into radiative channels (xx ~^ 

77 or XX -»■ 

we wish to extract 

how production/destruction of nuclei Hi due to photo-dissociation depends on DM mass. As in the 
case of hadro-dissociation, the amount of photo-dissociation is proportional to the amount of DM 
annihilation, dn x /dt. The amount of photo-dissociation also depends on the number of photons 
produced in an annihilation event and the subsequent electromagnetic cascade. Because there is 
no mass gap for the photon, the number of photons produced is proportional to the visible energy 
(e^ and 7) from the annihilation event, E v [ s [33J. For example, annihilation to electrons or muons 
have E v { s = 2m x and E v i s ~ 2m x /3, respectively. 

The amount of photo-dissociation implies that the Boltzmann term for production/destruction 
of nuclei Hi due to photo-dissociation is inversely proportional to the DM mass, 

dn Hi 
dt 

Therefore, if DM annihilates into radiative channels, constraints on (av) from BBN scale as m x , 
which is confirmed by precise numerical calculations for m x > 10 GeV [16[ fl"8| [T9] and shown in 
figure [TJ This scaling can be expected to hold down close to the binding energy of 4 He, E b A-R e ~ 20 
MeV. In the case of non-thermal production of 6 Li this scaling holds down to ~ 60 MeV. This is 
because the daughter 3 H or 3 He from photo-dissociation of 4 He needs about 10 MeV kinetic energy 
to efficiently capture on background 4 He [8], thus requiring E^ > 60 MeV. Note that constraints 
from the CMB also exhibit the scaling behavior (av) oc m x |20l \21\ [22j [23"1 [24"t [25] since the amount 
of visible energy from the annihilation is what dictates how many photons will be produced. 

3 Constraints from BBN 

Relic dark matter annihilations alter the primordial abundances of nuclei from standard BBN 
predictions. In the case of s-wave annihilation, as discussed in the previous section, the scaling 
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Figure 1: BBN constraints on hadronic (blue) and radiative (green) annihilation channels from D and 
3 He overproduction, respectively. The solid blue (green) line is our estimate based on the scaling behavior 
(<jv) oc m i- 5 ± - 2 (((T-u) oc m x ) for hadronic (radiative) annihilation. The long dashed lines indicate bounds 
on XX ~* bb (blue) and XX ~^ e. + e~ (green) from Hisano et. al. |16) . The dotted blue and green lines 
is a rendering of where we expect these scaling laws to break down, but their exact shape and placement 
requires precise numerical treatment. The red line is the bound from the CMB, whose normalization is 
set by constraints on XX ~^ bb at 95% C.L. from Natarajan [25] and extended using the scaling (av) cx 
m x (20] [21] [22j [23[ [24] [25] . For CMB, the important quantity is visible energy injection, and therefore 
annihilation to quarks or charged leptons give roughly the same constraint. 



of these effects on DM mass and annihilation cross-section is given in equations ([5]) and ([6]) for 
hadronic and radiative annihilation channels, respectively. 

In this work, we use BBN to place an upper bound on (av) for the previously unconsidered low 
DM mass region of MeV < m x < 10 GeV. To do this we adopt the bound on (av) at m x = 100 GeV 
from Hisano et. al. [16] and extrapolate the bounds to the low mass region using ([5]) and ([6]). 

While a precise treatment requires numerical simulation, results from other works |16 |. I1T |. IT8] . I19j 
indicate that the scaling behavior in ([5]) and ([6]) is robust. We will discuss where this scaling 
behavior breaks down and how precise our estimates are. 

Since we adopt the constraints of Hisano et. al. [16] as a starting point, the observational 
abundances used to place bounds are the same as in that work. These observational abundances 
are discussed in appendix lAl 

3.1 Hadronic constraints 

As discussed in equation ([5]), constraints on (av) from hadronic energy injection are approximately 
proportional to fn x ~ a where a ~ 0.5 comes from the dependence of the nucleon multiplicity on 
y/s = 2m x ; (njv) is approximately proportional to m". As m x approaches a GeV, this scaling is 
expected to break down and constraints on (av) to rapidly weaken since annihilations no longer 
produce energetic nucleons that dissociate background 4 He. 

The scaling (av) oc m x 5 assumes the dominant effect on primordial abundances is hadro- 
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dissociation of background 4 He. While hadro-dissociation does dominate, there may be other 
processes which contribute a small but non-negligible amount. For example, DM annihilations 
may further enhance primordial D, 3 He, and 6 Li abundances through (1) photo-dissociation of 
background 4 He by energetic photons produced in the hadronic jet and from electromagnetic energy 
loss and (2) more background 4 He to dissociate from p — n interconversion by pions at T > 100 
keV. While these contributions are sub-dominant to hadro-dissociation of background 4 He, precise 
numerical treatment is needed to evaluate the relative size and scaling with m x of all possible 
effects. The scaling behavior with m x might also be modified by the number of secondary energetic 
nucleons produced from hadro-dissociation processes, although this effect is subdominate [9]. 

Our goal in this work is to derive conservative estimates. Therefore, to parametrize our igno- 
rance on the interplay of all effects in calculating constraints on (av ) for a given mass, we include 
an uncertainty in the scaling 

{av) oc m^' 5±5 (7) 

where we take the uncertainty to be 5 = 0.2 [16]. 

In Ref [16] the authors consider annihilation into b quarks, XX ~~ ^ bb. From the primordial 
deuterium abundance they derive an upper bound (av) < 7.7 x 10 -25 cm 3 /s for m x = 100 GeV. 
As shown in figure [TJ the scaling ([7]) implies (av) < (av)th = 3 x 10~ 26 cm 3 /s for 

m b <m x < 11.4+! 3 , GeV. (8) 

For annihilation to bottom quarks, the constraint is cutoff at mj, ~ 5 GeV. While the above 
constraint is for XX ~ > bb, similar results are anticipated for annihilation to other quark species. 
However, in the case of annihilation to lighter quarks, we expect the constraint on (av) falls off as 
m x — > mjv ~ GeV, but where and how rapidly this occurs can only be obtained from numerical 
analysis. 



3.2 Electromagnetic constraints 

Constraints on (av) from electromagnetic energy injection are proportional to m x . Electromag- 
netic energy primarily alters primordial abundance through photo-dissociation of 4 He, leading to 
enhanced production of other nuclides. Therefore, the scaling (av) oc m x holds as long as the DM 
annihilation and subsequent shower produce photons energetic enough to dissociate 4 He. 

As discussed earlier, the observable consequences of overproduction as a result of He photo- 
dissociation are largest for 3 He. In |16| . Hisano et. al. study the annihilation channel XX e + e~ 
and derive an upper bound of (av) < 3.9 x 10~ 23 cm 3 /s for m x = 100 GeV from 3 He overproduction. 
As shown in figure [TJ the scaling (av) oc m x implies (av) < (av) th = 3 x 10" 26 cm 3 /s for 

30 MeV < m x < 77 MeV (9) 

The constraints fall off when the electrons from the DM annihilation can no longer produce 
photons energetic enough to dissociate 4 He. The photo-dissociation cross-section for He, a^^^ , 
proceeds through the giant dipole resonance and only becomes efficient when EL > 25 MeV [8[ 134] . 
An electron from the DM annihilation event with E e ± = m x near this threshold produces energetic 
photons via inverse-Compton scattering off background photons. Inverse-Compton scattering is a 
relatively hard event, with ~ 80% of the momentum transferred to the scattered photon. Therefore 
we expect the constraints on XX ~^ e + e~ from 3 He overproduction to fall off around 30 MeV. 
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3.3 Other annihilation channels 



The constraints ([SD and ([U]) for hadronic and electromagnetic energy injection are obtained from 
DM annihilation to bb and e + e~ , respectively. It is useful to consider how these constraints would 
change for different annihilation channels. Nucleosynthesis is sensitive to the amount of hadronic 
and electromagnetic energy injected into the bath and is therefore essentially decoupled from the 
high energy details of the dark matter annihilation. Therefore, in order to understand how different 
DM annihilation channels alter BBN, it is sufficient to specify the annihilation products and ask 
what particles are produced in the resulting hadronic and/or electromagnetic cascade. 

In the case of annihilation into hadronic channels, we do not anticipate the constraints from 
different quark species (u, d, s, c) to change much from the case of annihilation to bottom quarks. 
Injected hadronic energy primarily alters BBN through dissociation of nuclei by energetic nucleons 
resultant from the hadronization of the initial state quarks. At energies much larger than the 
nucleon mass, the spectrum of nucleons does not significantly depend on the initial state quark. 
As the DM mass gets near the nucleon mass, E > m^, the nucleon spectrum from hadronization 
will differ between initial state quark species. For lighter quark species, the main change from the 
constraint for XX bb §5§ is that the bounds will fall off at a DM mass lower than ~ 5 GeV. 
Therefore, for annihilation to quarks, we generically expect (av) xx ^, q q < (av)th for few GeV < 
m x < 11 ± few GeV where the lower limit comes from the inability to produce energetic nucleons 
in the annihilation event and will have a slight dependence on initial quark species, beyond the 
precision of these estimates. 

Annihilation to photons, electrons, or muons will alter BBN through the injection of electro- 
magnetic energy. From equation ([6]) , constraints on the DM annihilation cross-section are inversely 
proportional to the amount of injected visible energy (electrons and photons), 

(av) oc m x f^-J . (10) 

In the case of annihilation to photons or electrons E v [ s = 2m x . For annihilation to muons , -Evis rs "' 
2m x /3 and therefore the constraints on {av) for XX ~ > A t+ / U_ are weakened by a factor of ~ 3 from 
the case of XX ~^ e + e~ . Therefore, we estimate that 3 He overproduction bounds (w) xx _^^+ „- < 
1 x 10 _25 (m x /m^) cm 3 /s and is valid for m x > m^. 

BBN constraints on annihilation to taus, XX ~~ ^ t + t~, come from the injection of electro- 
magnetic energy from the r decay products. Despite the large hadronic branching fraction in 
the r decay, since m T < 2m at, annihilation to r's cannot produce nucleons which would alter 
nucleosynthesis through hadro-dissociation of 4 He. While the hadronic products of r decay (pi- 
ons and kaons) will slightly affect the primordial 4 He mass fraction through interconversion of 
background n and p at 100 keV < T < MeV, they predominantly alter standard BBN through 
the electromagnetic energy in their decay products. Taking the average fraction of visible en- 
ergy in a r decay to be 0.31 |16| . overproduction of 3 He from photo-dissociation of 4 He bounds 
(av) xx ^ T+T - < 1 x 10- 24 (m x /10 GeV) cm 3 /s. 

Finally, we comment on annihilation to neutrinos. Injection of energetic neutrinos during nu- 
cleosynthesis has been considered in the context of a long-lived particle decay [35], [36] [37] . While 
neutrinos have impact nucleosynthesis much less than colored or charged particles, here we esti- 
mate the magnitude of effects on nucleosynthesis from annihilation to neutrinos and show that the 
effects are, at most, as large as 0(1O -4 ) the effect of radiative annihilation channels. High energy 
neutrinos may pair produce charged leptons off of background neutrinos, v + Pbkg ~~ ► l + + l~ • If 
this reaction happens after T ~ 0.5 keV, the photons in the electromagnetic shower of the charged 
leptons can photo-dissociate nuclei and alter their primordial abundances [35] [36] EZ] ■ Charged 
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leptons are also generically produced in the annihilation through radiation of a (real or virtual) 
weak boson by one of the final state neutrinos, e.g. XX ~ * VihW^ — > v^iVjlj^ 

Both processes — high energy neutrinos pair producing off the neutrino background and anni- 
hilation to a multi-body final state — are important for charged lepton production, and therefore 
for BBN constraints. For charged lepton pair production off a background neutrino the rate of 
interaction is approximately F vj}b _>i+i- ~ Gpm x T v nu,bkg- This reaction can occur as long as the 
initial neutrino is above the threshold for pair production, E v E Uhk = m x T u > mf, which for pro- 
duction of e + e~ occurring after T ~ 0.5 keV requires m x > GeV. The fraction of neutrinos that 
pair produce after T ~ 0.5 keV is then approximateljll 



/ 

Jt( 



Note that if the initial neutrino is energetic enough, high energy secondary neutrinos from neutrino- 
neutrino scattering can play a role [66\ 137]. As for the multi-body final state, the relative amount 
of charged lepton production via neutrino radiation of a weak boson is approximately G 2 F m x / '16ir 2 
for m x <C Mw and can be as large as a w /Air ~ 3 x 10~ 3 when the weak boson is radiated 
on shell. Thus, depending on the energy of the initial neutrino, both the high energy neutrinos 
themselves and branching to multi-body final states can have the same order-of-magnitude effect on 
nucleosynthesis, which is what was found in [37J. Based on these considerations, we estimate that 
bounds on annihilation to neutrinos are weaker than bounds on radiative annihilation channels by 
O(10~ 4 — 10 -6 ) for m x >GeV and by ~ G 2 F m x /lQir 2 for m x < GeV. For this reason, annihilation 
to neutrinos is essentially negligible for BBN. 



4 Conclusions and Future Work 

In this work we have examined the effects of dark matter annihilations during the epoch of big 
bang nucleosynthesis. We emphasize that, in terms of the annihilation itself, the magnitude of 
these effects depend only on the rate of energy injection {i.e. the DM annihilation rate) and the 
type of energy injected (hadrons, charged leptons, etc.). With this procedure, we have explained 
how changes to primordial abundances of nuclei scale with the dark matter mass and annihilation 
cross-section. These scaling behaviors are robust and have been observed in precise numerical 
treatments of DM annihilation during BBN [T6 | W7 \ \TE [ [T9] . 

The dependence of changes to nucleosynthesis on the dark matter mass and annihilation cross- 
section, along with results from precise numerical calculation [16], have allowed us to estimate 
constraints on (av) for the previously unconsidered low DM mass region MeV < m x < 10 GeV. 
Interestingly, our estimates indicate that BBN rules out generic s-wave annihilation to quarks 
(radiative e+e",77) for few GeV < m x < 10 GeV (30 MeV < m x < 500 MeV). 

Our results have focused on the case that the thermally averaged annihilation cross-section is 
independent of time (s-wave annihilation). For scenarios in which (av) depends on time, changes to 
the time-independent case can be understood from a modification of the rate of annihilation, r ann = 
n x (t)(av)(t). For example, if the annihilation is s-wave suppressed, (av) will decrease in time [26J 

2 Of course, if the primary annihilation channel is through neutrinos, annihilation to l + l~ is generated at loop 
order and is model-dependent. Given a model, this annihilation can be calculated and the bounds for 
used. We focus on radiation of a weak boson giving a three- or four-body final state since it is model-independent. 

3 Cosmological redshift for the energetic neutrino can actually occur since T vi>bk ^^i+i- < H. However, to not 
complicate our point, we do not include it in the estimate. Cosmological redshift lowers the fraction of neutrinos that 
pair produce charged leptons, but does not change our conclusions. 
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and therefore have less impact on BBN. On the other side, as studied in [19] , Sommerfeld |38| or 
Breit-Wigner |39j type enhancements lead to (av) increasing in time and having a stronger impact 
on nucleosynthesis. 

Besides making our estimates precise, our results warrant a full numerical treatment, which is 
currently in progress [3D], for a few reasons. Our estimates indicate that s-wave annihilation to 
quarks is ruled out starting around m x ~ 10 GeV, which is competitive with the current bound 
from CMB [23JE5]. In estimating constraints for light DM, we extrapolated the results of Hisano 
et. al. [16], whose philosophy was to provide conservative constraints. Given that BBN provides 
an independent, and possibly more stringent, constraint than CMB, it is worthwhile to perform a 
full numerical calculation that also quantifies the confidence limit on the constraints placed. 
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A Observational abundances of light elements 

Since we use the results of Ref [16] as a starting point for our analysis, we adopt the observed 
primordial abundances used in their analysis. 

The primordial D/H abundance is inferred from the QSO absorption line in metal poor systems, 

(n D /n H ) P = (2.82 ± 0.26) x 10" 5 . (12) 

This value, used in Ref [16j, is the weighted average of six observed QSO absorption systems. 
A more recent work by the same group includes a seventh measurement that doesn't change the 
central value, but lowers the dispersion to ±0.20 x 10~ 5 [19]. 

An upper bound on the primordial 3 He abundance is obtained from 3 He/D measurements in 
protosolar clouds, 

(n3 He /n D )p < 0.83 + 0.27. (13) 

While we do not use Li to place bounds in this work, we comment here on its measurement and 
the so-called lithium problem. 7 Li is observed in the atmospheres of metal-poor Population II stars 
in our galactic halo. Recently, 6 Li has been observed in the these systems as well [41] . although 
the total number of systems with observable 6 Li is somewhat controversial due to observational 
difficulties in distinguishing 6 Li and 7 Li spectra [42J. 

Both the 7 Li and 6 Li measurements conflict with the theoretical values predicted by standard 
BBN and are collectively referred to as the lithium problem (for a recent review, see [43] ). For Li, 
the observed abundance is lower by a factor of about three than the theoretical value predicted 
by standard BBN, with a 4 — 5<t significance. In the case of 6 Li, the observed abundance is more 
than a factor of 10 2 larger than the standard BBN prediction. While the solution to the lithium 
problems may be unknown, a conservative approach, as used in [16] , is to take the observed 6 Li and 
7 Li abundances as upper limits for their primordial abundances. There is also an extra uncertainty 
in the primordial abundances that comes from the possibility that stellar burning may deplete 
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primordial 6 Li and 7 Li [33]. In this case, the conservative approach is to include this depletion as 
an uncertainty that raises the primordial abundance |16j . 
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